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Abstract 

We reanalyze the effect of magnetic fields in BBN, incorporating several fea- 
tures which were omitted in previous analyses. We find that the effects of 
coherent magnetic fields on the weak interaction rates and the electron ther- 
modynamic functions (p e , P e , and dp e jdT^) are unimportant in comparison 
to the contribution of the magnetic field energy density in BBN. In conse- 
quence the effect of including magnetic fields in BBN is well approximated 
numerically by treating the additional energy density as effective neutrino 
number. A conservative upper bound on the primordial magnetic field, pa- 
rameterized as ( = 2eB rms /(T^), is ( < 2 < 0.27p u ). This bound can 
be stronger than the conventional bound coming from the Faraday rotation 
measures of distant quasars if the cosmological magnetic field is generated 
by a causal mechanism. 
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1 Introduction 



It has frequently been suggested that a significant cosmological magnetic 
field is present 0, ||, U and over the last few years several proposals have 
been made how to generate such a field @]. Big bang nucleosynthesis (BBN) 
provides us with a tool that can be used to probe physics in the early universe 
and as such is invaluable in the study of primordial magnetic fields. In 
fact, BBN appears to be the only tool presently available that can explore 
magnetic fields in the pre-recombination universe and hence the physics that 
goes into generating the magnetic field. 

The effect of magnetic fields on BBN was discussed in the early pioneering 
works of Greenstein || and Matese and O'Connell ||. The conclusions and 
understanding in these papers were broadly correct. But now, with advances 
in our understanding of BBN and faster computers, we can make tighter 
quantitative estimates of the light element abundances ( 4 He, D and 3 He). 
The issue has also been examined in the recent literature, first by Cheng, 
Schramm and Truran (CST) [0] and subsequently by Grasso and Rubinstein 
(GR) ||. Our analysis differs from these calculations in that we have included 
several factors omitted in these works, corrected some apparent errors and 
resolved some of the conflicts. We also provide analytic arguments based on 
the Euler-Maclaurin expansion to deduce the changes in various quantities 
due to the magnetic field. 

There are two primary effects of magnetic fields on BBN : (i) the magnetic 
field energy density contributes to the cosmological expansion rate, and, (ii) 
the electron phase space changes due to the magnetic field. A priori, we 
cannot say if the dominant effect that determines any changes to the light 
element abundances is (i) or (ii), so it is necessary to go through the full 
analysis and to evaluate the magnitudes of both effects. The determination 
of the effects of the electron phase space change are quite involved and it is 
here that we disagree with CST and GR. However, after the full evaluation 
is carried out, the dominant effect of the magnetic field turns out to be its 
direct contribution to the cosmological expansion rate. The computationally 
expensive contributions to the weak rates and the electron thermodynamic 
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functions can be neglected, and in practice the numerical BBN computation 
can be avoided altogether by expressing the bound on the magnetic field 
energy density as an equivalent bound on the effective number of neutrinos, 
a number commonly derived in the BBN literature. This should be contrasted 
with the claims of GR who find that the change in p e dominates, and CST 
who find that the changes in the weak interaction rates dominate. We agree 
with the original conclusions of Matese and O'Connell. 

Although a rigorous comparison of the effects (i) and (ii) can only be 
made after a systematic evaluation, one way to understand why the effect 
on the expansion rate is more important than the change in the interaction 
rates is to note that the expansion rate changes in proportion to B 2 while the 
interaction rates change in proportion to aB 2 , with a ~ 1/137 being the fine 
structure constant. Any other numerical factors occuring in the interaction 
rates turn out not to be large enough to overcome the suppression by a. 

Our analysis assumes a magnetic field that scales as: 

BocR- 2 = bT 2 (1) 

where R is the scale factor of the universe and T v the neutrino temperature. 
Our analysis also assumes that the cosmological expansion is well described 
by the Friedmann-Robertson- Walker (FRW) model. The latter assumption is 
only justified if the magnetic field is sufficiently tangled on scales smaller than 
the horizon - otherwise we would need to consider the anisotropic expansion 
of the universe. At the same time, the assumption in (H) is valid only on 
scales large enough for the plasma conductivity to keep the magnetic field 
frozen-in. (Today this frozen-in scale is of order 10 13 cms.) Recent work 



by Jedamzik et. al. |L5[ considers the magnetohydrodynamical evolution of 
the system and concludes that the magnetic field can dissipate on yet larger 
scales which would correspond to scales of several megaparsecs today. The 
BBN constraints that we derive on primordial magnetic fields are therefore 
only valid for magnetic fields that are coherent on scales larger than the scale 
on which magnetic fields can dissipate and smaller than the horizon scale. In 
particular, magnetic fields generated during inflation are likely to be coherent 
on super-horizon scales and have not been fully dealt with by our analysis. 
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As the dominant effect of magnetic fields on BBN is due to the contri- 
bution of the field to the energy density, it will be convenient to express the 
energy density in magnetic fields in terms of the number of equivalent ad- 
ditional neutrino species. The energy density contributed by the magnetic 
field is 

R 2 

where, B 2 is the volume average of \B(x)\ 2 . The contribution of N v light 
(m <C lMeV) neutrino species to the energy density of the universe is 

7vr 2 

Pu = —N u Tt. (3) 
p 120 " y ' 

and so ps may be written as a number of equivalent additional neutrino 

species 
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where b = B/T 2 is the constant which we seek to constrain. This gives 



^ = (4) 



Pu + Pb = {K + AN?)p v {N„ = 1). (5) 

For standard BBN N v = 3. Krauss and Kernan (KK) argue that a very 
conservative bound on AN U is AN V < .8 0, leading to the bound b < 3.4. 

The bound in KK on N u is derived by requiring that at least 5% of 
1000 Monte Carlo runs simultaneously satisfy the abundance constraints 
D+ 3 He/H< 10~ 4 and Y v < .25, for any value of r], the baryon to photon 
ratio, for a given number of neutrinos. The random variables in the the 
BBN Monte Carlo uses 11 nuclear reaction rates and the neutron lifetime as 
random variables. The result, N u < 3.8, is conservative in that the abun- 
dance limits assumed allow for generous systematic errors in the observed 
4 He/H abundance in HII regions, and the likelihood of anomalous presolar 
abundances of D and 3 He. 

In the following sections we lay the physics groundwork for our bound 
on primordial magnetic fields. First we consider the effects of the B field on 
the electron phase space, then the resultant changes in the weak interaction 
rates and the time-temperature relation. Finally we tie these all together 
and numerically solve for the neutron fraction as a function of temperature. 
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2 Electron Phase Space 



The dispersion relation for the electron in the presence of a uniform magnetic 
field is 

E = (p 2 e + m 2 e + eB(2n + s + l)) 1/2 (6) 

with e the electron charge, B the magnetic field strength, and n the prin- 
cipal, and s the magnetic quantum number of the Landau level. This re- 
lation is valid for a magnetic field that is small compared to the critical 
field B c = m 2 e /e and has only been derived at zero temperature. For stronger 
fields and at high temperatures we expect further corrections to this formula. 
For example, the next order correction would be to include the anomalous 
magnetic moment of the electron. However, the analyses available in the lit- 
erature |10| indicate that these higher order corrections generally have a very 



small effect at the magnetic field strengths and temperatures of interest and 
so we will only consider the dispersion relation in @ throughout this paper. 
The full treatment of the problem is computationally impractical until an 
analytic approximation is developed for E(n, s,p, B). 

After summing over the electron spin in the zero-chemical potential 
phase space of electrons appearing in weak interaction rates and electron 
thermodynamic functions changes according to 

^f F o(E ) ^ l { 2-S M ) d ^f FD (E B ,T) (7) 

with E\ = p 2 + m 2 + 2eBn and fFD(E,T) the Fermi-Dirac distribution at 
temperature T: 

We have neglected the chemical potential, <ft of the electrons since 0/T 7 is 
of order the baryon to photon ratio, rj, which is small( O(10 -10 )); hence any 
B and ^-dependent correction to the standard result would be second order 
in small quantities. The electron thermodynamic functions are given by the 
usual prescription of integrating the relevant quantity over phase space. Thus 
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eB 



n,, 



+°° r+OO 

E(2-M/ f F D(E B ,T,)dp z (9) 

-DO J -°° 



(2vr) 2 f 



(2vr; 



E(2-5no)/ E B f FD (E B ,T^)d Pz (10) 



-DO " "°° 

e £ +2? „ . r+°°El-ml 



^ e = 7f^E(2-^o)/ -^^fFD(E B ,T 7 )dp z (11) 

(12) 



n e does not appear in the time-temperature relation, but is useful for illus- 
trative purposes, e.g. the perhaps counter-intuitive result that the electron 
number density increases as the magnetic field increases. This is easy to show 
by expanding @ in an Euler-MacLaurin expansion. 



T 7 u r p 2d p 



n e ~ -U4 



4vr 2L 7o ! + eV /p 2 +^l/T2 

+ (Cr, 7 ) 2 r dry 1 (13) 

24 7o Jtf + m 2/ T 2 ! + cosh ^2 + m 2/ T 2 

+ ^((Cr, 7 ) 4 )] 

where £ = 2eB/T^ = 2eb, and rv 7 = (T^/T 7 ) 2 . Although we will quote limits 
on 6 or (, which are constants, we will also use z = C r wy, which is what 
usually appears in the equations. Clearly from flUf), 

lim ^ > (14) 
b->o dB 

The interpretation is that there are two effects of increasing the magnetic field 
strength on the electron phase space - 1) the energy cost of populating n > 
levels becomes greater due to the larger spacing of the Landau levels, and, 
2) the areal element in the phase space decreases due to the decreased cross- 
sectional area of each Landau level. The second effect is stronger than the 
first and causes the number density of electrons to increase with B. For p e , 
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the sign of the effect is less clearcut (indeed it depends on the temperature): 



■[4 



2* r ?yJi?+mllT*d P 



4tt 2L Jo i + e ^P 2 +rnl/T* 

+ 1 _ Vy^Tg7g\ 

12 Jo \2{1 + cosh ^jrf + ml/T*) i + e Vp 2+m */ T * J 
+ 0(z% 

The first term inside the square brackets gives the electron energy density 
when B = 0. The second term is the correction to the energy density due to 
the presence of a magnetic field and can be negative or positive, depending 
on the temperature T 7 . 

Assuming that the integral in the z 2 term of fll5]) is 0(1), it is clear that 
the change in the energy density of electrons in the field is smaller than the 
pure magnetic field contribution to the energy density. Since the former is 
due to coupling of the magnetic field to the electrons, it is proportional to 
z 2 oc aB 2 times small numerical factors of order 0.01 — 0.1. This is smaller 
by 0(a) than the B 2 /8tc term. We expect the phase space corrections to the 
weak interaction rates to also be of order z 2 times small numerical factors. 

In p, only the changes in p e were taken into account; but the changes in 
P e and dpe/dT^ must also be taken into account since they appear in the BBN 
time-temperature relation. We will spare the reader the explicit formula for 
dpzjdT^ which is quite a lengthy expression with many extra terms involving 
dB/dT^. 



3 Neutron Fraction 

The binding energy per nucleon of helium-4, the most strongly bound light 
nucleus, is 7.1 MeV. Yet, well below this temperature the abundance of 4 He 
compared to hydrogen is quite small, because the high energy tail of the 
photon distribution breaks up any deuteron (Q = 2.2MeV) which forms. At 
a (slightly flsh^ dependent ) temperature of T 7 ~ .07 MeV the high energy 
tail of the photon distribution is no longer energetic enough to break up all 
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the deuterons. When this "deuterium bottleneck" is passed, all the available 
neutrons are swiftly incorporated into 4 He, with trace amounts of other light 
nuclei, D, 3 He and 7 Li produced. To a very good approximation, Y p , the 
primordial 4 He mass fraction is given by twice the neutron fraction at the 
breaking of the deuterium bottleneck. Thus to discover the importance of the 
effects of the corrections to the BBN dynamics, for a first pass it is sufficient 
to examine the evolution of the neutron fraction with all strong interations 
turned off. 

With respect to comparison with observed abundances the 4 He produc- 
tion is the most likely to be noticeably affected by the changes in BBN 
dynamics in the presence of the magnetic field. This is because the 4 He ob- 
servations are sensitive at the 5% level or better, whereas all the other obser- 
vations are sensitive only at the 50-100% level. Also, changes induced in the 
nuclear reaction rates when they are important (T 7 < 0.07MeV) are likely to 
be much smaller than those induced in the weak interaction rates when they 
are important (T 7 ^> 0.07MeV), both because the electron magnetic moment 
is much larger than the nuclear magnetic moment, and because the magnetic 
field decreases as ~ T^, and is therefore small by the time the nuclear 
reactions take place. Finally, the other elements are sensitive to changes 
in the time-temperature relation due to the magnetic field only because of 
the very small change in the temperature at which their number density 
falls out of quasi-equilibrium, whereas the neutron fraction at the breaking 
of the deuterium bottleneck (and hence the 4 He abundance) is sensitive to 
the change in the equilibrium neutron fraction during weak equilibrium, the 
change in the departure of the actual neutron fraction from it equilibrium 
value, the change in the freezeout temperature of the weak interactions, the 
change in the fraction of neutrons which decay between weak freezeout and 
the breaking of the bottleneck, and the change in the temperature at which 
the deuterium bottleneck breaks. 

The neutron fraction evolves according to 

dX 

^ — — A rap X n + \ pn (l — X n ) (16) 
where \ np is the rate per nucleon of neutron conversions to protons and \ pn 
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is the rate of proton conversions to neutrons. 



with a corresponding formula for X pn . For the numerical calculation we left 
out the 3-body reaction, pe~v — > n which never becomes important. If 
the weak interactions were in equilibrium, which would be the case if the 
interaction rates were much faster than the Hubble expansion rate, then the 
neutron fraction would assume its equilibrium value 

X e g = (lg) 
^np ' Ap n 

The weak interaction rates at these low energies are given simply by 
integrals of p e E e E 2 over the available phase space. For example, 

Kr^ru = ° fT ^ V ^ 39A) z E(2-<U r dp z E 2 u g(E e /T,)g(E u /T u ) (19) 

Here Gf is the Fermi decay constant, gy and qa are the vector and axial- 
vector coupling constants. The electron energy is given by 

E e = (p 2 e + m 2 e + 2eBn) 1/2 , (20) 

the neutrino energy by 

E v = A - E e , (21) 

with A = 1.293MeV being the neutron-proton mass difference. The sum in 
( |T9"D is taken to a maximum of n c =[(A 2 — m 2 e )/2eB\ (where [x] is the largest 
integer < x) , and the momentum integral is up to p m = A 2 — ml — 2neB. 
g is the Fermi blocking factor and is given by: 

g (E/T) = 1 - f FD (E/T) = (eM-E/T) + l)" 1 . (22) 

The other weak interaction rates have a similar form, with appropriate factors 
of fpD and g in the integrands, and with the limits on the sum and integral 
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given by the requirement of positive energy for the electron and neutrino. 
For example, 



A 



G 2 F T 2 (g 2 v + 3g 2 A ) 
(2tt) 3 



n=0 



dp z E 2 g(E e /T^)f(E u /T„) (23) 



With the exception of beta decay we find that the magnetic field effects 
on the electron phase space decreases all the weak rates. For small z we can 
see this once again using the Euler-McLaurin expansion: 



A 



G 2 F T 2 (g 2 v + 3g 2 A ) 



nu^pe 



(2tt)* 



dx 



z 2 T 2 
12 





' 1 dG-(E e )' 


'Pmin(0) 


_E e dE e 



Pmin (■*■) 

dp z 



G-(E e (p z ,x))dp z (24) 



x=0 



where 



P 



E e (p z ,x) = (p 2 z + m\ + x) 1/2 
(A 2 



2 rr2\l/2 

m e - xT 7 ) ' 



and 



Thus 



G-(E e 

1 dG_(E e ' 
E P dE P 



(E e - A) 1 



Ep - A 



1 



1 _|_ e (E e -A)/T u I + e -B e /T 7 



1 



1 



E, 



1 _|_ e (E e -A)/T v I _|_ e -B c /T 7 



(25) 
(26) 

(27) 
(28) 



E„- A 



1 



T 



+ 



Ep 



l + e ~(Ee-A)/T v Tl + e E e/ T ' 



We see that only the second term in the expression in square brackets is 
negative and that (for x = 0), by the time E e is large enough for the term 
to dominate over the other two terms in the square brackets, the prefactor 
is exponentially suppressed. The coefficient of the z 2 term of (|24|) is there- 
fore negative, and the B-field decreases this rate. Similar arguments can be 
applied to the other 2^2 rates. Though the conclusions are not always so 
clear analytically, our numerical studies show a decrease in all the 2-body 
rates for B > 0. 
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For beta decay the momentum cutoff introduces an n cutoff, and with 
B > B c = (A 2 — m 2 e )/2e only the n = term survives in the sum. Thus it 
is easily seen that \ n ^, pe - v increases linearly with B for B > B c . We find 
numerically that it also increases (though not linearly) for B < B c . 

There are three competing effects for the weak interactions. At high tem- 
peratures the weak rates are slower. This by itself does not imply anything 
about the equilibrium neutron fraction, since the neutron fraction (see fll8f )) 
depends on a ratio of rates; but, because X pn decreases more than X np , there 
is a small decrease in the equilibrium neutron fraction (see Figure 1). This 
is the first effect. However, below about 1 MeV, the weak interaction rates 
are too slow to maintain this equilibrium. Since X^ q is decreasing with tem- 
perature, the actual neutron fraction lags behind, and is therefore always 
greater than, the equilibrium value. Because the B field lowers the interac- 
tions rates, X n is less able to track X^ q , and so lags behind the dropping 
equilibrium neutron fraction more than usual. This second effect tends to 
increase the neutron fraction at a given temperature. Finally at later times, 
T 7 <,.2MeV, neutron decay dominates the weak interaction rates. The B 
field causes the neutron decay rate to increase, lowering the neutron frac- 
tion. However, since neutron decay is the dominant reaction only at lower 
temperatures, the B field has already decreased significantly, and the change 
in the decay rate is negligible. 

If these changes in the reaction rates were the only effects, then the decline 
in the equilibrium fraction and the increase in the neutron decay rate at low 
temperatures would outweigh the increased departure from equilibrium (due 
to the lower reaction rates) at high temperatures. The net result would be 
a lower neutron fraction and hence a lower 4 He abundance as can be seen 
from the dotted line in Figure 2. 

However, as we shall discuss in detail below, these are not the dominant 
effects. 



11 



4 Time- Temperature Relation 



To solve for BBN abundances we must follow not just the evolution of the 
neutron fraction, but also the evolution of the photon and neutrino temper- 
atures with time. Initially, T 7 > 5MeV, the neutrinos and electromagnetic 
fluids are assumed to be in perfect thermal contact and Tl l = T % v . Later, the 
photon and neutrino temperatures evolve differently due to electron-positron 



annihilation |TT, [T2|], which heats the photons, but not the neutrinos. This 



assumes zero net energy exchange between the neutrino and electromagnetic 



fluids, a very good approximation for BBN |fL3|j . Assuming that the neutri- 
nos are decoupled well before this annihilation takes place, T u evolves as the 
inverse scale factor 

RT V = constant. (29) 
For the photons we use energy conservation in the electromagnetic plasma. 

d[( Pem + Pem)R 3 } = R 3 dP em , (30) 

with p em = p e + p 7 and P em = P 7 + P e . It follows that 

dpem/dT^ = -3(p em + Pe^R^dR/dT^ (31) 

The magnetic field energy density is not to be included in (^). This would 
be appropriate only if the large scale B field was in thermal equilibrium with 
the electromagnetic fluid, and hence evolved according to T 7 rather than the 
inverse scale factor. 

The Hubble expansion law for the radiation dominated early universe is 



R dt V 3 y J 

with pT representing all forms of energy density in the universe. (Note that 
the magnetic field energy is included here but not in (|3~1~D). From (|3l|) and 



(|32"D we obtain the time-temperature relation which we are seeking, 



dT'y 3H^ em ^ >ern f33") 

dt dpem/dT^ 
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The analogous equation for the neutrino temperature, from fl29|) and (^) is 

§ = < 34 > 

The complication in solving (|3^) in the case of a magnetic field is that p e 
depends both on T 7 and on B = bT 2 , p e = p e (T y , T V ;Q. We therefore write 

dp e dp e dp e dT u 



dT^ dT 1 dTy dT y 



(35) 



If we now define 



§MMW |f = ^X) (36) 
we can rewrite (|33|) Using ( |3~4l) and (|35|) as, 

dt dp^/dT^ + h 3(p em + P em ) ' ' ' ' 



The right-hand sides of equations (0) and (|37|) depend only on T U1 T 7 , and 
the constant £ (or equivalently 6) and can be integrated together to obtain 
T„(t) and T 7 (t). 

Although the details must be obtained numerically, the basic picture is 
clear. The increased energy density due to the B 2 /8n term means that at any 
given temperature (T 7 or T„) the universe is younger and expanding faster 
than in the standard scenario (B = 0). This has two important consequences. 
First, the weak interaction rates, which are already slower than in the B = 
case have an even harder time maintaining the equilibrium neutron fraction. 
Since the equilibrium neutron fraction is dropping, this increased lag will 
cause the actual neutron fraction to be higher than for B = 0. Second, the 
time between weak interaction freezeout and the breaking of the deuterium 
bottleneck is shortened. Thus there is less time for the extra neutrons to 
decay than if B were zero. This also acts to increase the final neutron fraction. 
Although there is a lowering of the temperature at which the bottleneck is 
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broken, this is a second order effect and the change in this temperature 
turns out to be only about 0.1% and is therefore ignored. 

The net result is that if one takes into account only the change in the 
time-temperature relation, the final neutron fraction rises and so does the 
4 He abundance. 



5 Results 

We are left to determine whether the neutron fraction rises because of the 
change in the time-temperature relation, or falls because of the changes in 
the weak interaction rates. We find that the effects of the change in the 
time-temperature relation from the magnetic field energy density are much 
larger than those in the time-temperature relation or the weak interaction 
rates due to the changes in the electron phase space. We have not been 
able to demonstrate this convincingly using only analytic methods; however, 
basically it is because the changes in the phase space are proportional to 
aB 2 , whereas the change in the energy density is proportional to simply to 
B 2 (with no a). 

We integrated the coupled differential equations (0), ( ft7| ) and ( |TE| ) to 
obtain the photon and neutrino temperatures and the neutron fraction as a 
function of time. We ran the computer code from T 7 = 5MeV to .07MeV, 
when the deuterium bottleneck is broken and almost all the remaining neu- 
trons are converted to 4 He. We did this for several values of £ = 2eB/T 2 . 

In Figure 2 we display the evolution of the neutron fractions for ( = 0, 
and for £ = 2 in each of five different limits. The long dashed curve shows 
the ( = case; the solid curve shows X n with the full calculation; in the 
dotted curve we have turned on only the dependence of the weak interaction 
rates on the magnetic fields; in the short-dash curve we have only turned 
on the dependence on the magnetic fields only of the electron thermody- 
namic functions; and with the dash-dotted curve we have only turned on 
the dependence of the expansion rate on the energy density of the magnetic 
field. It is clear that the latter is by far the dominant effect. The change 
in X n (T 7 = .07MeV) from the canonical case, B = 0, is —2.0 x 10~ 5 for 
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weak interactions only, < 5.0 x 10~ 6 for electron thermodynamics only and 
+4.8 x 10~ 3 for magnetic field energy density only. Thus the limit on B is 
equivalent to one on the number of neutrinos, and for N v < 3.8, our limit on 
B is given by, 

C < 2 (38) 
In units of the neutrino energy density this can also be expressed as 

Pb < .27 Pu , (39) 

where p v is the energy density contributed by the standard 3 light (<< 1 
MeV) neutrinos. 

An additional point of interest is that the change in the electron ther- 
modynamics during the era of e + e~ annihilation causes a slight decrease 
in T u /T y with respect to the canonical value, (4/11) 1 / 3 . This can be un- 
derstood heuristically as follows. The magnetic field increases the electron 
number density. This is equivalent to some fraction of an additional degree of 
freedom residing in the electromagnetic fluid in the early universe. When the 
e + e~ pairs annihilate the additional entropy is transferred to the photons, 
thus "heating" them slightly more with respect to the neutrinos. For the 
values of the magnetic field allowed by BBN the change is small enough that 
it has little effect on other cosmological parameters which depend on the neu- 
trino temperature, such as the mass of a light neutrino required for Q u = 1. 
If very large magnetic fields were allowed, for example, with ( = 8.0 (16.0), 
we find 

= -1.7% (-5.6%), (40) 

and this effect becomes significant. 

In Figure 3 we display the evolution of the neutron fraction as a function 
of T 7 from T 7 = 3 to T 7 = 0.07 for ( = .1,2,4. The points on the curve 
where t = .1, 1, 10 and 100 seconds are marked. The two dominant effects 
which we described are manifest. At high temperature, the 2 2 weak 
interactions are dominant. The higher the B field the faster the universe 
expands at a given T 7 , and the more the actual neutron fraction lags behind 
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the decreasing equilibrium abundance. At low temperatures, neutron decay 
is dominant. The higher the B field, the sooner the universe cools sufficiently 
to break the deuterium bottleneck, so the sooner neutron decay terminates 
and the remaining neutrons are converted almost instantaneously to 4 He. 
This is seen by the earlier termination of the X n line. Both these effects lead 
to a higher neutron fraction. 



6 Discussion 



Our analysis of the effects of magnetic fields on BBN have allowed us to place 
constraints on the field strength present at the BBN epoch. The analysis is 
based on the following assumptions: 

a) The cosmological model is an FRW model. This means that the mag- 
netic field does not lead to anisotropic expansion of the universe. For this 
the magnetic field should not be coherent on the scale of the horizon such as 
happens in proposals that generate magnetic fields from inflation. However, 
it seems likely that if the magnetic fields are not too strong, the constraint 
we have derived using an isotropic FRW universe should still apply. If there 
are significant departures from FRW universes, the magnetic field would be 
constrained by measurements of the isotropy of the universe stemming from 
measurements of the microwave background [117 . 



b) The magnetic field is assumed to have a scaling given by (0) and hence 
we have assumed that it does not dissipate. The shortest scale that does 



not dissipate has been argued to be given by the magnetic Jeans length [15 



Another scale that is probably relevant is the scale on which the magnetic 
fields are frozen in due to the high electrical conductivity of the cosmological 
plasma Jl|, || . 

c) The dispersion relation is assumed to be given by (|]) which is known 
to be incorrect when the magnetic field is strong (compared to m 2 e /e) or at 
high temperatures. But current analyses [1(J indicate that the corrections 
are small in the range of temperatures and magnetic field strengths that are 
relevant for us. We assume that these corrections will not change our result 
significantly. This assumption can be removed once we have a better analytic 
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grasp of the corrections to the dispersion relation as it would be computa- 
tionally infeasible to use the high magnetic field strength and temperature 
dispersion relation that is known in its present form. 

Within the framework of these assumptions, we have derived a BBN con- 
straint on the present cosmological volume averaged magnetic field strength 
Brms ( equivalent ly, on the magnetic field energy density) defined by 

BLs = ^J v d 3 x\B(x)\ 2 . 

Using flSgp and the evolution equation ([!]), the constraint is: 

B rms < 10- e G . (41) 

We would like to emphasize that this constraint is a local constraint and not 
a constraint on the magnetic field strength at any particular coherence scale. 

An oft quoted constraint on primordial magnetic fields comes from the 
Faraday rotation measures of distant quasars. This constraint is 

B < 1(T 9 G (42) 

where B is the line-averaged magnetic field: 

B = - [ ds- B(x) 
L Jy 

and the curve V is the line of sight (null geodesic) of length L from the quasar 
to the earth. 

To compare ( f4"T| ) with (f4"2"|), we need to know something about the spec- 
trum of the magnetic field |16| . If the field is homogenous on the quasar 
distance scale, then a direct comparison can be made and we see that the 
BBN constraint is much weaker. This comparison is relevant if the magnetic 
fields arise from inflation where it is assumed to be coherent on horizon scales 
(subject to the discussion in point (a) above). If, on the other hand, a causal 
mechanism is found to generate a magnetic field, the field is likely to be co- 
herent on scales that are much smaller than the horizon, say on a scale £. 
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For example, in the case of magnetic fields generated in cosmological phase 



transitions, we expect 16 



J-'rms — 




and the Faraday rotation constraint in (013) translates into 



B rms < y —10 9 G 

The BBN constraint (|41|) can therefore be much stronger than the Faraday 
rotation constraint if L/£ is large. 

Finally we would like to project the results of this paper from a different 
viewpoint. This is that the current astrophysical constraints on magnetic 
fields constrain the large scale magnetic field strength. That is, they say 
something about the Fourier modes with wavelengths comparable to quasar 
distances. BBN, on the other hand, provides a complementary tool for con- 
straining the magnetic field energy density which includes contributions from 
all wavelengths. And for causal mechanisms that generate magnetic fields, 
BBN is likely to probe the magnetic field spectrum at the smallest wave- 
length while the Faraday rotation measure is mostly sensitive to wavelengths 
of order the quasar distance scale. 
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Figure Captions 



Figure 1: The equilibrium neutron fraction X^ q as a function of the photon 
temperature T 7 for several choices of the magnetic field evolution parameter 

c 

Figure 2: For ( = 2, the neutron fraction X n near the nucleosynthesis temper- 
ature region, O(70) keV, is shown for separate pieces of the full calculation. 
The changes from the electron thermodynamic functions, (p e , P e , dpe/dT^ ), 
the weak interaction rates X n ^ p , and the magnetic field energy density p#, 
on the neutron fraction are shown individually. Also plotted is the neutron 
fraction with all/none of the above corrections. 

Figure 3: The neutron fraction X n as a function of the photon temperature 
T 7 for several choices of the magnetic field evolution parameter (. 
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